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PRIMORDIAL BLACK HOLES AS DARK MATTER: ALMOST ALL OR ALMOST NOTHING 
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ABSTRACT 

Primordial black holes (PBHs) are expected to accrete particle dark matter around them to form ultracompact 
minihalos (UCMHs), if the PBHs themselves are not most of the dark matter. We show that if most dark matter 
is a thermal relic, then the inner regions of UCMHs around PBHs are highly luminous sources of annihilation 
products. Flux constraints on gamma rays and neutrinos set strong abundance limits, improving previous 
limits by orders of magnitude. Assuming enough particle dark matter exists to form UCMHs, we find that 
^pbh ^ 10 -4 (for /«dmc 2 ~ 100 GeV) for a vast range in PBH mass. We briefly discuss the uncertainties on our 
limits, including those due to the evolution of the UCMH luminosity as it annihilates. 

Subject headings: dark matter — early universe — diffuse radiation — gamma rays: diffuse background 



1. INTRODUCTION 

The early Universe was extremely smooth, but it is possible 
that there were rare but large perturbations. Any large per- 
turbations (Sp /p > 0.3) collap sed to form Primordial Black 
Holes (PBHs; Hawking 1971). PBHs are expected to ac- 
quire a halo from the surrounding particle dark matter back- 
ground, if_fhej^_themselTCS do_jiot make up most of the dark 
matter dMack et alj 120071: iRicottil 120071: IRicotti eTal] 120081: 
iRicotti & Gouldll2009|). These halos are called Ultracompact 
Minihalos (UCMHs: IScott & Sivertssonll2009t) . 

The abundance of PBHs with masses above 1000 M is 
strongly constrained (SIpbh ^ 10~ 8 ), using the ef fects of bary- 
onic accretion on the CMB energy spectrum (IRicotti et alJ 
120081) . and that of masses below ~ 10 15 M Q is constrained 
by the n on-observation of gamma rays from PBH evapora- 
tion (see ICarr et alj 12010, and references therein). In the 
planetary-stellar mass regime, there are modest constraints 
(^pbh 5. IP" 1 ) from microl ensing ( e.g.. lAlcock et alj[l 998. 
I200U iTisserand etaLll2007L but see lAlcocket all 120001) and 
the dynamic s of structure fo rmation and widely separated bi - 
nary stars dYoo etalj|2004t but see also lOuinn et all 12009). 
Over a huge range of possible PBH masses (10~ 15 - 10~ 9 M©), 
the c onstraints in their abundances are weak or non -existent 
(e.g.. lSeto & Cooravll2007T:lAbramowicz et alj|2009h . 

The dark matter UCMHs around PBHs allow new ways 
to con strain the abundance of PBHs. Ricotti & Gould 
(120091) suggested searching for the microlensing signature 
of UCMHs, also pointing out that there could be a gamma- 
ray signal f rom dark matter annihil ation in the dense core of 
the UCMH. IScott & Sivertssonl (120091) developed this further, 
considering the gamma-ray signal from individual, nearby 
UCMHs of various masses, and their detectability with Fermi. 

We develop strong new constraints on dark matter annihila- 
tion in UCMHs around PBHs, using general assumptions that 
can be evaded only in exotic dark matter models. These limits 
apply if PBHs do not make up all of the dark matter, so that 
enough particle dark matter exists to form UCMHs. We strive 
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for order-of-magnitude accuracy. 

2. DARK MATTER ANNIHILATION 

We assume dark matter is a thermal relic, following anni- 
hilation freezeout in the early Universe, consisting of Weakly 
Interacting Massive Particles (WIMPs). We assume it may 
self-annihilate in the late Universe (i.e., there is no large 
WIMP-antiWIMP asymmetry), and that the final states are 
Standard Model particles. For thermal relic dark mat- 
ter, r^DM^ 2 = 3 x 10~ 27 cm 3 s _1 /(oAv), relating the dark mat- 
ter d ensity to its thermall y-averaged annihilation cross sec- 
tion dKolb & Turner! I1990I) . We take ft D M = 0.3 and h = 
ff /(100 km s" 1 Mpc -1 ) = 0.7. We consider masses above 
the GeV range, probed at colliders, and below a PeV, set by 
the unitarity bound. We use the relic (oav), assuming it inde- 
pendent of v; this holds for standard s-wave annihilation. Our 
results cannot be evaded by reducing (o>v), as this would lead 
to an unacceptably large finM- 

While (ctav) is fixed, the particular Standard Model final 
states are not; however, strong con straints on dark mat ter an- 
nihilation can be set in any case dBeacom et al.1 120071) . We 
then limit the abundance of UCMHs from constraints on dark 
matter annihilation, and thus the abundance of PBHs them- 
selves. 

Constraints on dark matter annihilation from gamma-ray 
signals are well known. Direct annihila tion produces a 
spectral line at wdmc 2 (Mack et al.1 [2008). Annihilation 
into quarks, charged leptons, or gauge bosons produces a 
gamma-ray continuum when those particles decay. Electron- 
positron final states produce gamma rays, through inter- 
nal bremsstrahlung, with branching fraction Br(7) » a w 
0.01. This process is an electromagnetic radiative correc- 
tion; it is independent of surrounding matter density, oc- 
curs in any process involving charged particles, and pro- 
duces gammaj^s_wWi_energies up to m DM c 2 dBeacom et al.l 
l2005t iBell & Jac ques 2009). Additional radiation is produced 
through synchrotron and Inverse Compton energy-loss pro- 
cesses. 

It might be thought that neutrinos are invisible annihila- 
tion products; in fa c t, there are stringe nt bounds on them 
dBeacom et all 120071: lYuksel et al.1 120071) . The sheer size of 
neutrino telescopes like IceCube compensates for the small 
detection cross section, and the atmospheric neutrino back- 
ground falls off steeply with energy. Other Standard Model 
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annihilation products are intermediate to the above cases. 

3. UCMH PROFILES AND LUMINOSITIES 

According to iMack et all (120071) and iRicottil (120071) . a 



UCMH (with or without a PBH) at z \ 
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UCMHs with PBHs at matter-radiation equality have M eq = 
Mpbh- The dark matter mass density within R u -(z) is the same 
no matter what the PBH mass is, with a number density n n - ra 

where mmn = otdmc 2 /(100 GeV). 



1.8 x 10 5 cm 3 at z, 



eq> 



Within 7? eq = ^n-feqX the UCMH density profile goes as p oc 
r~ 3 / 2 , because the PBH dominates the mass; at larger radii, a 
r~ 9 / 4 density profile holds dBertsch inger 1985). This density 
profile is shallower than in lRicotti & GouldT (l2009). which as- 
sumed p oc r -9 / 4 throughout, and that adiabatic contraction 
increased the dark matter density further. 

WIMPs have nearly radial orbits around the PBH, sampling 
a wide range of densities in each orbit. Thus, the dark matter 
annihilating in the inner regions of the UCMH is often stored 
much farther out. To answer whether WIMPs with apocenters 
of R eq survive to the present day, we calculate the number of 
annihilations a WIMP is expected to experience during the 
Universe's history. This is the number of annihilations per 
orbit times the number of orbital periods in a Hubble time: 

(AW= 2 (y)/ n DM (r)(a A v)^de, (2) 

where we have defined 6 = as apocenter and P is the orbital 
period. A WIMP can only annihilate once, so the density pro- 
file is valid at present only if (N ann ) < 1. Assuming the orbit 
is Keplerian and nearly radial, 
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where r a is the apocenter radius and v a is the WIMP's tangen- 
tial velocity at apocenterQ Using Eq. [TJto find R eq and 
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from the bulk velocity dispersion in iRicotti & Gouldl (2009), 
we find 



(AU,) » 0.03 m 



100 



lOGyr 



1-0.005 In 



PBH 



Mr. 



(5) 



Therefore, WIMPs with apocenters of R eq mostly survive to 
the present, despite venturing much further in. 
The accretion and resultant WIMP orbits at R eq must be 

nearly radial: ctdm *C \J '2GMpbh / 'Req- Using v a (Eq. |4]l for 
the initial velocity dispersion ctdm, we find accretion is ra- 
dial for all reasonable PBH masses. WIMPs also have resid- 
ual thermal dispersion, v t h e rm ra 1.3 cm s _1 m7o (z/zeq) _1 ; con- 
servatively, we find radial accretion at R eq holds for Mpbh > 



4 We approximated K(k) = 
ln(7TV C F/2)], where k 2 = 2e/(e 
if k 2 -C —1, which holds for nearly radial orbits. 



jQ /2 d9/Vl-k 2 sm 2 e as l/v^Fp + 
1), and e is the eccentricity. This is valid 



2 x lO~ 15 mjQ O (z/z e q)~ 3 M0. Accretion at later times than z ra 
3000, at smaller radius than R &q , or from the low end of the 
WIMP velocity distribution might still form a UCMH; so our 
limits may roughly hold for smaller PBHs. Finally, WIMPs 
must miss the central PBH. WIMPs with an apocenter r a have 
apericenter r p ra v 2 r 2 /(2GMpBH); if the apocenter is R eq , then 
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Thus, WIMPs with apocenters of R eq remain outside the 
Schwarzschild radius if M PB h < 1740 M Q . 
The UCMH annihilation luminosity is L ann = 



f ci ' 2tt r 2 n(r) 2 ( a Av)moMC 2 dr, or 

T 9 M 2 BH c 2 (a A v) 



32-7T HlDM^e, 



In 
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where r m i n is a minimum-radius cutoff. Even if WIMPs with 
small apocenters annihilate away by the present, WIMPs on 
radial infall orbits with larger apocenters generate a time- 
averaged density profile of p oc r" 3 / 2 between apocenter and 
pericenterQ Therefore, we use r p eq as r m ; n . The UCMH lumi- 
nosity is 



24 Lr. 
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ignoring a very small logarithmic term in Mpbh- The linear 
scaling with Mpbh arises because the UCMH mass scales lin- 
early with PBH mass, and the WIMP density (and annihila- 
tion lifetime) at R eq is always the same. When p oc r" 3 ' 2 , each 
decade in minihalo radius produces the same annihilation lu- 
minosity; the total luminosity therefore is not strongly depen- 
dent on the WIMP orbits' eccentricities. The eccentricities 
themselves depend weakly on M PB h, since r a eq / r p eq oc Mpg H . 
Thus the total luminosity is understood as the mass of the 
UCMH halo (proportional to Mpbh) annihilating over a fixed 
annihilation timescale (set by n tT ) multiplied by some slowly 
varying logarithmic factor accounting for the UCMH's inner 
regions. 

A steeper density profile inside R eq would increase the lu- 
minosity but shorten the WIMP survival time. Eq. [2] im- 
plies, for WIMPs with a r~ 9//4 density profile and r a = R eq , 
(N, nn ) = 15OOm-i o (MpB H /M ) 08 , so t, nn ra 7 m - l 00 Myr. For a 
r 2 profile, (AU) = 22m7i o (MpBH/M ) 05 (z ann ra 10 when 
mioo = 1). Efficient annihilation introduces ~ SIpbh worth 
of high energy radiation into the Universe. At high red- 
shift, the CMB energy spectrum and reionization history 
severely con strain the gamma-ray and cosmic ray injec- 
tion r a te (e.g.JFixsen et alJl996t |Padmanabhan & Finkbeiner 
120051: Mapelli et alj I2006f) . At low enough redshift, the 
gamma-ray and neutrino backgrounds still constrain the an- 
nihilation products directly. Thus, it would be difficult to 
have SlpBH ra fJwiMP, unless perhaps the annihilation were en- 
tirely into neutrinos. Furthermore, such efficient annihilation 
would mean that present UCMHs around PBHs have depleted 
inner halos, with consequences for microlensing searches 

5 Consider a shell of non-interacting WIMPs at apocenter r a and with iden- 
tical tangential speeds. The time-averaged mass at each radius is proportional 
to the time they spend there: (dM(r)) oc dr/v r , where v r is the radial velocity. 
For nearly radial orbits, v r ss ^JIGMjr, except near pericenter and apocen- 
ter. Thus (dM(r)) oc r l l 2 dr\ since (dM(r)) = 4irr 2 (p(r))dr, (p(r)} oc r~ 3 / 2 
except near pericenter and apocenter. 
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dRicotti & Gouldl 120091) . More work on the effects of anni- 
hilation in steep density profiles is needed. 

3.1. Assumptions and Other Considerations 

This annihilation luminosity depends on some simple con- 
siderations. We assume that PBHs have not evaporated yet, 
which is valid for Mpbh ^ 10~ 18 M Q ; and that most of the 
dark matter is made of WIMPs and not PBHs, otherwise there 
would be too few WIMPs to form a UCMH. 

We require that the density profile along a WIMP's orbit not 
evolve significantly. Most importantly, we require that dark 
matter accretion after z eq does not increase the mass within 
R eq . Such accretion would initially increase the UCMH lu- 
minosity, but would also shorten the time WIMPs survive 
in the UCMH. We also assume that annihilation of WIMPs 
with small apocenters does not affect our calculations. As 
a limiting case, if WIMPs annihilate into a constant density 
core with p = p(R eq ), then L ann w O.42L (Mp B H/M )m^ o ; 
this weakens our bounds by a factor of ~ 60, but still im- 
plies fipBH <C 0.1 for otdmc 2 < 10 TeV. Furthermore, since 
the PBH dominates the mass within R eq , adiabatic contraction 
should be unimportant. 

Perhaps the most important other effects we do not consider 
are those of accreted baryonic material, which may cool radia- 
tively and collapse efficiently. Adiabatic contraction will not 
be important if the PBH dominates the mass within R eq . How- 
ever, baryonic matter can be optically thick to gamma rays, 
reducing the apparent UCMH luminosity. We estimate the op- 
tical depth as r ~ n e ajR t0i . If the mass of baryons within R eq 

is f b M DM « f h M PBH , then t rs O.S/^Mpbh/Mq) 1 / 3 . Klein- 
Nishina effects will reduce this optical depth. Thus for 
smaller PBHs (M < 10 M Q ), we can ignore baryonic opac- 
ity. Neutrino limits are unaffected by opacity. 

4. COSMIC GAMMA-RAY BACKGROUND CONSTRAINTS 

Annihilation of dark matter can produce gamma rays 
with significant power near todmc 2 - Gamma rays below ~ 
100 GeV contribute directly to the extragalactic background. 
Gamma rays with energy > 100 GeV for z > 1 cascade down 
in energy by pair production and Inverse Compton processes 
with ambient photons, contributing to the background at lower 
energies. 

The gamma-ray emissivity of the Universe at energies 
above 100 MeV is limited by EGRET observations to Q ms . x = 
8 x 10" 35 erg s" 1 cm" 3 dCoppi & Aharonianl 119971) . The 
number density of PBHs is then limited to be «pbh ^ 
<2max/ (£annBr(7)), where Br(7) is the branching fraction into 
gamma rays. If dark matter annihilates into charged parti- 
cles, there must be internal bremsstrahlung, with branching 
ratio Br(7) m a w 0.01, which we take as a minimum branch- 
ing fraction. This number density is converted into a limit on 
Upbh by multiplying by the mass of the PBH and dividing 
by the critical density, p c = 9.20 hjo x 10~ 30 g cm" 3 , so that 

^PBH < 2maxMpBH/l>AmnBr(7)): 



PBH < 1-9x10 



- 5 /Br( 7 )V 
V 0.01 J 



(9) 



The upper limits on SIpbh are essentially independent of PBH 
mass: the number of PBHs for a given 57pbh scales as Mp BH , 
but the luminosity of each scales as Mpbh. 

5. MILKY WAY GAMMA-RAY CONSTRAINTS 



B0 

o 



lllll ' >" 






6 - 








, i ... i ... i ... i ... i .. . 



10 §10 t m DM ° £ 



4 

(GeV)] 



FIG . 1 . — Upper bounds on the abundances of PBHs as a function of WIMP 
mass. Bounds on annihilation into gamma rays (black; Br(7) = 1) and elec- 
trons (grey; Br(7) = 0.01) are shown, as well as neutrinos (Br(i^) = 1) (blue). 
Cosmic background limits are solid and Galactic limits are dashed. Gamma- 
rays are the easiest final state to detect, while neutrinos are the hardest, and 
other Standard Model final states would give intermediate limits. 

Stronger constraints can be obtained by taking advantage of 
the higher-than-average number density of PBHs in the Milky 
Way. PBHs in the Milky Way are close enough that annihi- 
lation gamma rays are not attenuated; such gamma rays do 
not have to compete with the entire gamma-ray background 
above 100 MeV, but only with that near mdmc 2 - 

Suppose the density of PBHs, «pbhCs), tracks the dark 
matter density. Then the integrated gamma-ray inten- 
sity on a line of sight out of the Milky Way is 7 = 
4^rL an n(«PBH)Br(7) J S(s)ds, where 6 is the dark matter over- 
density over the average cosmic dark matter density and 
(«pbh) is the average PBH density in the Universe. Then the 
abundance of PBHs is limited by 



PBH 



< 



47r/ obs MpBH 



iannBr(7)p c / 5(s)ds ' 



(10) 



To find the background 7 b s the UCMH radiation competes 
with, we use the Fermi-m easured extragalactic background 
spectrum (lAbdo et al.ll2010h for E < 100 GeV. For annihila- 
tion into gamma rays or internal bremsstrahlung from charged 
particles, most of the power is within one log bin in energy of 
otdmc 2 dBell & J acques 2009). We therefore find 7 b s by inte- 
grating the gamma-ray background as such. 

We use an NFW density profile for the distribution of 
UCMHs in the Milky Way, 6(r) = ^(r/r.r'd + r/r s T 2 , with 
S s « 45000 and r s = 27 kpc dStoehr : et alj|2003l) . The line of 
sight integral is not sensitive to the distribution of UCMHs in 
the inner Galaxy. The integral is more like that for dark matter 
decay than diffuse annihilation, since the intensity is linearly 
proportional to «pbh('")- We consider a sightline aimed di- 
rectly away from the Galactic Center, where the uncertainty 
in the profile should have the least effect and the signal is 
smallest, for a conservative result. 

The dependence of SIpbh with WIMP mass is shown in 
Figure [T] for various final states. The Galactic gamma-ray 
bounds in Figure Q] (solid) include the cosmic background 
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bounds above 100 GeV. The annihilation luminosity falls as 
the WIMP mass increases (Eq.0. Since the competing extra- 
galactic background falls as E 2 dN/dE oc E~° 4 in the GeV to 
100 GeV range, the allowed SIpbh increases with WIMP mass. 
The Galactic signal constrains Opbh < 10~ 6 for a WIMP mass 
of 100 GeV even if Br(7) = 0.01, smaller than the cosmic 
bound (Eq.|9]). 

With these abundances, we can calculate the lower lim- 
its on the mean distance to the nearest PBH, Apbh = 
[3tt/(4(5(«pbh))] 1/3 , where the local 8 = 89000: 

A PBH >220pcm 10 / , (11) 

for WIMP masses greater than 100 GeV. This implies a 7-ray 
flux of 

<I> < 1.0 x 10 cm s m,J, 

7 ~ \ Mq J \ 0.01 J 

(12) 

6. NEUTRINO BACKGROUND CONSTRAINTS 

Unlike gamma rays, neutrinos do not cascade down in 
energy as they travel through the Universe, although they 
redshift. We integrate the atmospheric neutrino spectrum 
(or diffuse neutrino background limits above 100 TeV) 
from m DM c 2 /e to m DM c 2 dGaisser & Honda! [20021) . and re- 
quire the neutrino flux from UCMHs be l ess than this; 
other w ise, they would ha ve been detected dBeacom et al.1 
120071: lYiiksel et al.1 1 20071) The measured data is re - 
ported in lAshie et all d2005l). iGonzalez-Garcia et al.l (120061). 
lAchterberg et alj (120071). iHoshina et al.l (12008b . lAbbasi et all 
(2009). and lDe Young et al.l (12009). 

In Figure [1] we show the Galactic (solid) and cosmic 
(dashed) bounds on PBHs from neutrinos. The atmo- 
spheric neutrino background also steeply falls with energy 
(E 2 dN/dE oc E 3 ), so that the bound on fipBH is fairly 
energy-independent up to 100 TeV. 

7. CONCLUSIONS 

In Figure [2] we show that even with the uncertainties in 
our estimates, our constraints are powerful for PBHs inside 
UCMHs. PBHs of many masses are ruled out to very small 
abundances, for a wide range in WIMP masses. At high PBH 
mas ses (> 1 Mq), CM B constraints become more power- 
ful dRicotti et al.1 120081) : at very low masses, Hawking radi- 
ation limits are more powerful (e.g., ICarr et afl l2010l). Mi- 
crolensing cons traints are weaker , but make no a ssumptions 
about WIMPs dAlcock et alJfcOOltlTisserand et al.1120071) . For 
10~ 15 M Q <Mpbh i$ 10~ 9 M , ours are the only constraints on 
PBHs aside from 57dm- Our conclusion depends on the stan- 
dard assumption that most dark matter is a self-annihilating 
thermal relic. Our analysis d oes not apply if all of the dark 
matter is made of PBHs (e.g., Frampton 2009), because there 
will not be any WIMPs to annihilate. For the smallest PBH 
and WIMP masses (M PB h < lO~ 15 m7 ( 3 l() M ), radial accretion 
may not hold. PBHs could remain important fo r other rea- 
sons (e.g., as seeds of massive black holes; see Ma ck et all 
120071) which do not require them to have a substantial abun- 
dance. The combination of our results with previous limits 



imply that PBHs either make up almost all of the dark matter, 
or almost none of it. 

Our results could be improved by detailed simulations of 
how the UCMH evolves as its inner regions annihilate. An- 
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Fig. 2. — Upper bounds on the current abundance of PBHs of mass A/pbh 
for a WIMP mass of 100 GeV/c 2 (left of the transition in RgwrefTJ. The 
limits in this work (blue; solid where baryonic opacity and non-radial accre- 
tion can conservatively be ignored, and dotted if baryonic opacity is small 
and accretion is radial) are the most powerful limits over a vast mass range. 
See the mild caveats in the text. Also shown are microlensing limits (short- 
dashed-long-dashed), CMB limits from FIRAS (short-dashed) and WMAP3 
(long-dashed), and evaporation radiation limits (dash-dotted). See text for 
references. 

nihilation could lower the luminosity of the inner regions of 
the halo over time, weakening our bounds. Steeper density 
profiles than r" 3 / 2 may increase the annihilation luminosity, 
but shorten the lifetime of WIMPs so that the halo annihilates 
away before z = 0. We expect that some combination of con- 
straints on gamma-ray and neutrino backgrounds, reionization 
history, and the CMB energy spectrum will generally require 
a small SIpbh- 

Further improvements could come from more detailed stud- 
ies of the annihilation products and their detectability. When 
considering annihilation into charged particles, we considered 
only gamma rays from internal bremsstrahlung, but charged 
particles can themselves radiate, such as th rough Inverse 
Compton scattering (e.g., Cirelli & Panci 2009). We are con- 
fident the limits for WIMP masses ab ove 100 GeV can be im- 
proved by new observations. Finally. iRicotti & Gouldl (2009) 
suggest that some UCMHs may exist without PBHs, formed 
from weaker initial perturbations in the early Universe. Lim- 
its on dark matter annihilation in these UCMHs may strongly 
constrain these weaker perturbations. 

We thank A. Gould, J. Rich, M. Ricotti, K. Stanek, G. 
Steigman, and T. Thompson for helpful discussions. This 
work was supported by NSF CAREER Grant PHY-0547102 
to J.F.B. and an Alfred P. Sloan Fellowship to T. Thompson. 
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